Phases and phase transitions in a 17(1) x U(l) system with 9 = 2ir/3 mutual statistics 
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We study a (7(1) x (7(1) system with short-range interactions and mutual 8 = 27r/3 statistics in (2+1) 
dimensions. We are able to reformulate the model to eliminate the sign problem, and perform a Monte Carlo 
study. We find a phase diagram containing a phase with only small loops and two phases with one species of 
proliferated loop. We also find a phase where both species of loop condense, but without any gapless modes. 
Lastly, when the energy cost of loops becomes small we find a phase which is a condensate of bound states, 
each made up of three particles of one species and a vortex of the other. We define several exact reformulations 
of the model, which allow us to precisely describe each phase in terms of gapped excitations. We propose field- 
theoretic descriptions of the phases and phase transitions, which are particularly interesting on the "self-dual" 
line where both species have identical interactions. We also define irreducible responses useful for describing 
the phases. 



I. INTRODUCTION 

One of the hallmarks of topological quantum phases is that 
they have anyonic excitations, which can be viewed as parti- 
cles with statistical interactions. Examples include quasiparti- 
cles in the Fractional Quantum Hall Effect^, spinon and vison 
excitations in Z2 spin liquids*^— and excitations in a variety of 
interesting fractionalized systems^ - — It is also fruitful to ask 
about possible new phases that such particles can have, as a 
way to access proximate phases and phase transitions involv- 
ing topological quantum states J2r— 

Unfortunately, direct Monte Carlo studies are hampered by 
the sign problem. It turns of that some such systems allow re- 
formulations where they become free of the sign problem and 
can be studied using unbiased numerical approaches. Interest- 
ing questions include, for example, what phases can result if 
there are two species of particles with mutual statistics that are 
both trying to condense. In this work, we pursue such a study 
of the effects of a statistical interaction on a model of two 
species of integer- valued loops with short-range interactions. 
We are able to reformulate this model so that it can be stud- 
ied on a lattice using Monte Carlo techniques. Previously^ 
we studied a model with two species of loops and mutual it 
statistics, which is also of interest in effective field theories 
of frustrated antiferromagnets^2r~ and other areas £r— We 
would like to extend this to study systems with general sta- 
tistical angle 8. We have found that 9 = n is a special case, 
and the properties of such models are qualitatively different 
when 9 7^ ir. In this work we study 6 = 2tt/3, and the results 
should exhibit behavior similar to that for general 9 ^ tt. 

Our model can be precisely described by the following ac- 
tion: 
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The index r refers to sites on a cubic lattice (the "direct" 
lattice), and R refers to sites on another, inter-penetrating 
cubic lattice (the "dual" lattice)^— Ji^{r) is an integer- 
valued current on a link r, r + /t of the direct cubic lattice, 
J2n(R) is integer- valued current on a link R,R + ft of the 
dual cubic lattice. We use schematic vector notation so that 



FIG. 1 : The contribution to the partition function is multiplied by a 
phase e' e for each cross-linking of the two currents shown in the fig- 
ure on the left. If we change the direction of one of the currents and 
get the figure on the right, the phase is e~ lS . When considering sym- 
metries of our model for 9 7^ n, we must only consider operations 
that leave the relative orientation of the current loops unchanged. 



Ji and J 2 represent these conserved integer-valued currents, 
and V ■ Ji = 0, V • J2 = 0. In the partition sum, a given cur- 
rent configuration obtains a phase factor e 10 or e~ l6 for each 
cross-linking of the two loop systems, dependent on the rela- 
tive orientation of the current loops, as shown in Fig. Q] This 
is realized in the last term of Eq. (JT}, by including an auxil- 
iary "gauge field" p 2 , defined on the direct lattice, whose flux 
encodes the J2 currents, J2 = V x 02. 

Figure [2] shows the phase diagram for the model with 9 = 
2ir/3. When both t\ and t% are small we have a phase [labeled 
(0) in the figure] where there are only small loops. When t\ is 
large and is small, we get a phase [labeled (I) in the figure] 
where one species of loop has proliferated, while the other 
species has only small loops. Since our model is symmet- 
ric under interchange of t\ and t<x, we get similar behavior 
in phase (II). Since these phases do not have both species of 
loops occurring at the same time, the statistical interactions 
are unimportant. 

We now consider the region of the phase diagram where t\ 
and *2 are similar, in particular in this work we will often study 
the "self-dual" line where t\ =t-x. In this region, if we were to 
neglect the statistical interaction (9 = 0), we would have two 
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phases: a "gapped" phase at low t±, t 2 , where there are only 
small loops, and a "condensed" phase at high t\, t 2 , with pro- 
liferated loops in both the J\ and J 2 variables. The condensed 
phase would have two gapless modes, one from each species 
of loop. The transition from the gapped to condensed phase 
would be two decoupled XY transitions. If we turn on the 
statistical interaction, we find qualitatively different behavior. 
For small t\, t 2 we again get a gapped phase, but for larger 
ti, t 2 , we get a phase, labeled phase (IV) in Fig. [2] where 
the statistical interactions are manifest more dramatically. We 
will see below that in this phase both species of loop are con- 
densed, however there are no gapless modes. This phase is 
distinguished from phase (0) by a non-vanishing correlation 
between currents of different species. Such correlations were 
identically zero in the 9 = tt case, and this phase was not 
present in that model. 

If we increase t\ and t 2 still further, we get a phase, la- 
beled phase (III) in Fig. |2j which is a condensate of bound 
states composed of three particles in the J\ variables and an 
anti-vortex in the J 2 variables. This is a (2+l)-dimensional 
analogue to the 6*-term induced "dyon condensates" in (3+1) 
dimensions described in Refs. l32hl34l Loosely speaking, these 
composite states appear so that the system can avoid destruc- 
tive interferences from the statistical interaction. For exam- 
ple, the statistical interaction in Eq. (Q~|i is inoperative when 
the J-currents are present only in multiples of three, while the 
precise description of the phase (III) is obtained by employ- 
ing duality approaches in the main text. The transitions from 
phase (IV) to phases (0) and (III) occur at interesting multi- 
critical points, and we study the system's behavior at these 
points. 

The outline of this paper is as follows. In Section HH we 
reformulate the model, Eq. ([T), in a sign-free way so that 
we can study it in Monte Carlo. Section UTT1 contains the re- 
sults of the Monte Carlo study. These results are presented 
in terms of the correlation functions of the original J vari- 
ables of Eq. (Q]), which already allows us to distinguish all 
phases. In Section [IV] we introduce several additional exact 
reformulations of the model using duality transfor m 1 1 i 29 i 35 ^ 
summarized in the Appendix. These reformulations enable us 
to precisely describe each phase in terms of variables which 
are gapped in that phase This leads us to propose continuum 
field theories for the various phase transitions in our model in 
Section[V] In Section[VT]we derive equations for "irreducible 
responses" which provide a physical way to characterize the 
"condensates" that give phases (IV) and (III). We conclude in 
Section [VTIl bv comparing with the 8 = tt case and discussing 
further generalizations. 



n. MONTE CARLO METHOD AND MEASUREMENTS 

In Ref. [lj| we described a method of reformulating mod- 
els with short-range interactions and statistical terms, such as 
Eq. (HJ, in a sign-free way so that they can be studied in Monte 
Carlo. We review that method here, since in this work we have 
defined new measurements based on the sign-free reformula- 



tion. First, we pass from J\ variables to conjugate 27r-periodic 
phase variables by formally writing the constraint at each r: 



6[V- Ji(r) =0] = 



. exp[— i0 r (V • Ji)]. (2) 



To be precise in our system with periodic boundary condi- 
tions, we also require total currents of J\ and J 2 to vanish. In 
this case we can write J 2 = V x p 2 and the action (Q]i is inde- 
pendent of the gauge choice for p 2 . We enforce the zero total 
current in J\ with the help of fluctuating boundary conditions 
for the <j>-s across a single cut for each direction \i = x,y,z 
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This gives the following partition function: 
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where the action is given by: 

[V xp 2 (r)] 2 
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(5) 



- ^ VVillain [<Pr+n ~ 4>r ~ &P2^(r) - 7/A M ,o] 



Vvmain is the "Villain potential", which is obtained by sum- 
ming over the J\ variables: 



exp[-Vvillain(cMl)] 
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In the actual Monte Carlo, we use (f> r , 7^e(— tt, tt), 
p 2fl (r)eZ, and perform unrestricted Metropolis updates. One 
can show that physical properties measured in such a simula- 
tion are precisely as in the above finitely defined model. 

In this work, we monitor "internal energy per site", e = 
S/Vol, where Vol = L 3 is the volume of the system, and 
compute heat capacity, defined as 



C = (( £ 2 )-<6> 2 )xVol. 



(7) 



To determine the phase diagram, we monitor loop behavior 
by studying current-current correlations, which are defined as: 
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where a and b are the loop species and \i and v are directions; 
■Mfc) = Er Jap(r)e- ij; - p . We trivially have C*»(fc) = 
C°^(— k). Because of the vanishing total current, we define 
the correlators at the smallest non-zero fc; e.g., for we 
used k = (0,^,0) and ft = (0,0,^). In this work we 
are interested in the correlations between currents of the same 
species, (ft), also known as the "superfluid stiffness". For 
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example, C 22 can be measured easily in our Monte Carlo, 
since we have direct access to J2 = V x p2- 

We are also interested in the correlations between currents 
of different species, and we first need to find the correspond- 
ing expressions in terms of the variables in Eq. (|5}. We can 
couple the original J variables to external probe fields A ext 
by adding the following terms to our action: 



SS 



AT\r) +^J 2 (B). A c 2 xt (R). (9) ^ 



We carry the fields Af^ through the reformulation procedure 
and then take derivatives of the partition function with respect 
to them. We obtain the following expression for the correla- 
tion between currents of different species: 
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In the above equation, it is important to note that J\ and J2 are 
defined on different lattices. In order to work with them on the 
same footing in fc-space, it is convenient to define R = r' + d, 
where r 1 is on the direct lattice and d is the offset between the 
two lattices. We can choose any convention for this offset, and 
we choose d = (1/2, 1/2, 1/2), which means that the sites of 
the dual lattice are located at the centers of the cubes forming 
the direct lattice, and we use such "physical" coordinates for 
all sites when defining the Fourier transforms. For a given 
variable W(r) on the lattice whose sites are labeled by indices 
r, the quantity V x W(r) is defined on a dual lattice. Now that 
we have defined the relation between the two lattices, we can 
precisely define the meaning of the curl operation in fc-space, 

[VxW] p (k) = 2ie pvfl sm(K/2)e ik " /2 e- ik ^ 2 W fl (k). (11) 

We can use symmetry arguments to determine some prop- 
erties of the correlators C"„(k). For simplicity, in this work 

we define k to be in the z-direction, so that fc = (0,0, k z ), 
and we only need to consider fj,, v <E {x,y}. For a sym- 
metry operation to leave our action in Eq. ([1) unchanged, it 
must preserve the relative orientation of two cross-linked cur- 
rents, like those in Fig.Q] One symmetry that satisfies this re- 
quirement involves mirror reflection about a plane while also 
reversing the direction of one species of loop. C^(fc) and 
C^ 2 (fc) change sign under such an operation about a plane 
perpendicular to the x-axis, and therefore must be zero. We 
can also use such an operation about a plane perpendicular to 
the z-axis to show that C^ 2 (fc) is an odd function of fc, and 
hence C 21 x (k) = C 12 y {-k) = -C} 2 (k). Our action is also 
unchanged if we take its complex conjugate while also revers- 
ing the direction of one species of loop. We can use this, along 
with our precise definition of the offset between the two lat- 
tices and of the Fourier transforms, to show that all the corre- 
lators (fc) are real. Lastly, we can use the ir/2 lattice ro- 
tation symmetry of the action to show that C""(fc) = (fc) 
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FIG. 2: The phase diagram for the model in Eq. ([T} with 6 — 2ir/3. 
Phase (0) contains no loops. Phase (I) contains proliferated loops 
in Ji and no loops in J2, while in phase (II) the variables are inter- 
changed. In phase (IV) both species of loops are condensed in single- 
strength loops. Phase (III) is a condensate of bound states comprised 
of three charges from one species of loop and a vortex from the other 
species. The precise description of these condensates is given in the 
text. 



and Cj 2 (fc) = —Cy 2 (k). Whenever we present numerical 
data, we have performed appropriate averages over all direc- 
tions to improve statistics. 



III. RESULTS 



A. Mapping out the phase diagram 



We determined the different phases of the model by looking 
at the stiffness C* 22 (fc) = C 22 (fc) = C 22 (fc), defined at fc = 
fcmin = (0, 0, 2n/L). Its L — > 00 limit is non-zero in phases 
(II) and (III) and vanishes in the other phases. Since our model 
is exactly symmetric around the self-dual line, we know that 
C 11 (fc m i, 1 ) is non-zero in phases (I) and (III). We found the 
locations of the phase transitions more accurately by studying 
C 22 (fc m i rl ) • L crossings. We took data in sweeps across the 
phase diagram (see Fig.|2]i, and defined the intersection of the 
C 22 (fcmin) • L curves to be the location of the phase transition. 
An example of such a sweep is shown in Fig. [3] The dots 
on the phase diagram in Fig. [2] are the locations of the phase 
transitions determined in this way. In all C 22 (k m - m )-L sweeps, 
we found that the crossings did not drift with increasing L, 
which suggests that these phase transitions are second-order. 

Let us now consider some limiting cases. The model with 
t\ = is a model containing only one species of loop. 44 
Our value for the position of the (0)-(II) XY transition (£2 ~ 
0.333...) is in agreement with prior work on this model4^ The 
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FIG. 3: A sample of the C 22 (fc m i n ) ■ L crossings from which the lo- 
cation of the phase transition was determined. Error bars come from 
comparing runs with different random seeds. Here the transition is 
from phase (IV) to phase (II). The data is taken for t± = 1.0, and we 
determine the phase transition to be at ti = 1.80 ± 0.01. There is 
no sign of drift in the location of the crossings, suggesting a second- 
order transition. Note that the J2 variables are condensed in both 
phases (IV) and (II). Nevertheless, C 22 (fc m i n ) • L ~ 1/L in phase 
(IV) and C 22 (fcmin) • L ~ L in phase (II), so the crossing analysis on 
the plot detects the transition. This subtlety is discussed in the text. 



transition is in the 3D XY universality class also for small, 
non-zero t\. 

For t\ — > 00, the Villain weight © vanishes except for 
a = 2irx (int), which enforces J 2 — V x a 2 = 3x(int). 
Therefore, at t\ — >• 00 the (I)-(III) transition is a transition 
from no loops of J 2 to J 2 loops of strength 3. One expects that 
this transition is XY-like, and similar to the (O)-(II) transition, 
but due to tripled J 2 , it should occur at a i 2 value nine times 
higher. We observed the (I)-(III) transition to occur at t 2 « 3 
for large t\, in agreement with this expectation. We give a 
precise description of phase (III) for finite t\, t 2 in Sec. [IV] 

In Fig. [4] we show C 22 (fc m i n ) • L along the self-dual line 
£1 = ^2, going through phases (0), (IV), and (III). Nei- 
ther of phases (0) and (IV) have a finite superfluid stiffness 
C 22 (fc m in), so to distinguish between them we use the corre- 
lator C;J 2 (fc m i n ), denoted as C 12 (fc m i n ) in what follows. A 
plot of C 12 (fc m i n ) along the self-dual line is shown in Fig. [5] 
C 12 (fc m in) ■ L vanishes in phase (0) in the L — >• 00 limit, but is 
non-zero in phase (IV), so the two phases are indeed different. 

We can understand the observations in phases (0) and (IV) 
as follows. The excitations in phase (0) are small loops in the 
J variables, which implies that in this phase C aa (k) ~ fc 2 
for small fc. For fc = fc m i n , this gives C 22 (fc m i n ) ~ 1/L 2 . 
The smallest excitation that contributes to C 12 consists of a 
small loop in each of the Ji and J 2 variables. An estimate of 
such contributions with cross-linking between the loops leads 
to C 12 (fc) - -sin(6>)fc 3 . 

In phase (IV), the J variables are condensed. One way of 
expressing this condensation is to replace the integer-valued 
J with real-valued variables j. This is equivalent to coarse- 
graining the model and integrating out the gapped vortices 



(see Sec. |V). If we define new gauge variables dji and dj 2 

V 2 ° j2 , then we can replace 
the original action Eq. (fTJ by an effective action in terms of 
the (in, oij2 variables: 



such that ji = V 2 ° jl and j 2 
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In the absence of the last "mutual Chern-Simons" (CS) term, 
this would be an action for two decoupled gauge fields, which 
would have two gapless modes. When the mutual CS term is 
included, it destroys the gapless modes. We can calculate the 
C aa (k) and C 12 (fc) correlators with respect to this gaussian 
action, and we find that C aa (k) ~ fc 2 ~ 1/L 2 for k = fc min , 
consistent with our data. We also find that C 12 (fc) ~ — i = 
— |5 = — j- for k = fcmin . This quantity is represented by the 
dotted line in Fig. [5] and we can see that our Monte Carlo data 
approach this value. 

Let us briefly remark on the use of C 22 (fc m i n ) • L crossings 
to determine the phase boundaries. It is natural to use these 
crossings on the (0)-(II) and (I)-(III) transitions, where we are 
going from a phase with only small J 2 loops to a phase with 
large J 2 currents. For the transition from phase (IV) to (II), we 
are going between two phases where the J 2 variables are con- 
densed. However, since in phase (IV) C 22 (fc m j n ) • L ~ 1/L 
while in phase (II) C 22 (fc m i n ) • L ~ L we can still use cross- 
ings in this quantity to determine the transition between the 
two phases. One might not expect, however, to see qualita- 
tively similar behavior between this transition and the transi- 
tions (0)-(II) and (I)-(III), yet this is what we observe. The 
reasons for this will be explained in Sec.QV] For the transition 
between phase (I) and phase (IV), C 22 (fc m i n ) ■ L ~ 1/L in 
both phases, and so we cannot use it to detect this transition. 



B. Transition (0)-(IV) along the self-dual line 

We now investigate the apparent multicritical points on the 
self-dual line. We first study the lower regime where phases 

(0) , (IV), (I) and (II) meet. We are interested in how the phases 
join. Due to the symmetry between t\ and i 2 , there are three 
scenarios, shown in Fig. [6] In a), all four phases meet at a 
point, while in b) there is a critical line segment on the self- 
dual line, and in c) such a segment is perpendicular to the self- 
dual line. Figure|7]shows C 22 (fc m i n )-i along the self-dual line 
near this transition. C 22 (fc m i n ) • L vanishes in phases (0) and 
(IV) [see also Fig.|H, but appears to have a finite value in the 
critical regime. If scenario c) were correct, we would expect 
C 22 (fc m in) ■ L to vanish at the (0)-(IV) transition since phases 

(1) and (II) should not influence its behavior. In addition, we 
have taken sweeps with t 2 = t\ + St, 5t = 0.002, which 
are lines parallel to the self-dual line and displaced from it by 
St. We found two distinct phase transitions near the critical 
point, so if scenario c) is accurate the line segment is < 0.004 
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FIG. 6: Different scenarios for how the phases can meet on the self- 
dual line. In scenario a), all four phases meet at a point. In b), phases 
(I) and (II) meet on a line segment on the self-dual line, and in c) 
phases (0) and (IV) meet on a line segment perpendicular to the 
self-dual line. We believe that scenario c) is unlikely, but we can- 
not distinguish between a) and b) at finite size. We can only state 
that if such a segment exists in scenario b), it is in the narrow range 
t = 0.335 - 0.35. 



FIG. 4: C 22 (fe m in) • L along the self-dual line. We can see that 
C 22 (fc m i n ) • L vanishes in phases (0) and (IV), despite its unusual 
behavior at the transition between these phases. C 22 (fc m i n ) ■ L di- 
verges in phase (III) due to the proliferation of triple-strength loops 
in the J variables. The vertical lines mark the locations of the phase 
transitions. More detailed data at the phase transitions is shown in 
Figs. |7]and[T0] 
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FIG. 5: C 12 (fc min ) • L along the self-dual line. C 12 (fc min ) • L van- 
ishes in phase (0), but it approaches a universal value in phase (IV). 
We predict this value to be at ~2tt/8 = —3, which is shown by 
a horizontal line. In phase (IV), C 12 (fc m m) • L approaches a non- 
universal value. The vertical lines mark the locations of the phase 
transitions. 



in size. For these reasons, we believe that scenario c) is not 
taking place. 

Furthermore, if there is a line segment as in scenario b), 
it is no larger than the region in Fig. [7] where C 22 ■ L is in- 
creasing with system size. We can further limit this segment 
by studying heat capacity shown in Fig. [8] and noting that the 
segment is no larger than the region where heat capacity in- 
creases with system size. We therefore estimate that the line 
segment is within the small range t G [0.335, 0.35]. Studying 
larger sizes could further narrow our bounds on the possible 
extent of the line segment, but at finite system size we cannot 



show that it does not exist. 

To determine the order of the transition at this point, we 
studied how the heat capacity increases with system size. We 
can see from Fig. [8] that the heat capacity has a sharp peak 
around t = 0.345. We also studied histograms of the total 
energy per site e. In the second-order case, these histograms 
would be singly-peaked, while in the first-order case we ex- 
pect to see two distinct peaks. An example of such a histogram 
analysis, taken at t = 0.345, L = 28 and 32, is shown in the 
inset for Fig. [8] We do not see two distinct peaks, however the 
histograms have a "flat top", suggestive of two peaks which 
are too close to be distinguishable on our finite sizes. This flat 
top suggests that we have a first-order transition. 

The data we have presented seems to suggest that we have 
a first-order transition in the form of scenario a), which would 
be highly unusual. We therefore propose two alternate scenar- 
ios. Firstly, the transition could be first-order and scenario b), 
with a line segment which is very small. However, in addi- 
tion to the small size of the line segment, if this scenario were 
true we would expect C 22 (k min ) to be finite at the transition, 
since it is non-zero in phase (II). Therefore we would expect 
C 22 (k m i n ) ■ L to increase linearly with L on the segment, and 
this is not consistent with Fig. [7] Alternatively, scenario a) 
could indeed be correct, and the transition could be second 
order. The behavior of C 22 (k m i n ) ■ L in Figs. [4] and [7] im- 
plies that we have very strong crossovers in our simulation 
variables: as we approach the transition from phase (IV), we 
need larger and larger sizes to see the eventual vanishing of 
C 22 (fcmin) • L ~ 1/iin this phase. It is possible that the un- 
usual shape of the energy histograms could be due to sampling 
in these variables. Studying the system at larger sizes could 
help to resolve these questions, by both more clearly resolving 
the histograms and further reducing the extent of the possible 
line segment of scenario b). 



C. Transition (IV)-(III) along the self-dual line 

Figure|9]shows the heat capacity in the regime where phases 
(IV), (III), (I), and (II) meet. The peaks in the heat capac- 
ity evolve only slowly with system size, suggesting a second- 
order phase transition. Figure [TOlshows the C 22 (fc m i n ) Lnear 
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FIG. 7: C 22 (fc m in) • L along the self-dual line at the bottom multi- 
critical point. C 22 (fcmin) • L vanishes as 1 /L for sufficiently large L 
in both phase (0) and phase (IV), but is non-zero at the transition. 
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FIG. 8: Heat capacity along the self dual line at the bottom multicrit- 
ical point. The sharpness of the peaks suggests either a small v or 
a first-order phase transition between phases (0) and (IV). The inset 
shows histograms of the energy per site e, at L = 28, 32, t = 0.345. 
The irregular shape of the histogram suggests that this is a first-order 
transition. 
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FIG. 9: Heat capacity along the self-dual line at the upper multicriti- 
cal point. The behavior of the peak with size suggests a second-order 
transition. The inset shows a histogram of the energy per site e, at 
L = 32, t — 2.66, which is where C has a maximum for this size. 
The single peak further suggests a second-order transition. We also 
studied histograms at all of the other points where heat-capacity was 
measured, and found only single-mode distributions everywhere. 
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FIG. 10: C 22 (fcmin) • L along the self-dual line at the upper multicrit- 
ical point. The crossings do not drift with increasing size, implying 
a second-order transition at i cr it = 2.62. 



this point. At this transition we are going from a phase with 
C 22 (fc m i n ) • L ~ 1/Lto C 22 (fcmin) • L ~ L, so we expect a 
crossing at the phase transition. We observe that this crossing 
does not drift with increasing L, further supporting the con- 
clusion that the transition is second order. Finite-size scaling 
arguments suggest that C 22 (fc m in) ■ L = f[(t — i C rit)£ 1/ ' I/ ] in 
our model. We can therefore try to collapse the C 22 (fc m in) ■ L 
data in Fig. [I0]to one curve by rescaling the horizontal axis 
by (i — £ cr it)L 1 / I/ . Applying this process, using £ cr ;t = 2.62 
inferred from Fig. [10] gives v = 0.8 ± 0.1, consistent with 
a second-order transition. We have also obtained histograms 
of total energy at all of the points in Fig. [9] and have found 
singly-peaked histograms at all points. The inset to Fig. [9] 
shows our data for t = 2.66, L = 32, which is the location 
of the heat capacity maximum in the figure. This phase tran- 



sition is a transition from a phase where the J variables are 
condensed in single strengths to a phase where they are pro- 
liferated only in triple strengths. However, we have used tech- 
niques for analyzing C 22 (fc m in) ■ L which are valid for the 
ordinary condensation of loop variables. This will be justi- 
fied by a more precise description of the two phases and the 
transition between them in Secs. HVI andfVl 



IV. ANALYSIS IN TERMS OF EXACT 
REFORMULATIONS 

Using the duality transform shown in the Appendix, we can 
derive several exact reformulations of the action in Eq. ([!} . We 
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can use these reformulations to describe each phase in terms 
of variables whose loops are gapped in that phase. The nature 
of the different phases and the transitions between them can 
then be described in terms of these variables. We can also 
introduce into our initial action external "probe" gauge fields 
coupled to J i and J 2 , by adding terms to the action similar to 
those in Eq. Q: 



SS 



k 



-k) ■ Al xt (fc) + < £V 2 . A| xt (k). (13) 



We can carry these gauge fields through the duality transforms 
as illustrated in the Appendix, and take derivatives with re- 
spect to them to obtain various exact relations between differ- 
ent current-current correlators. We will use such relations to 
better understand the behavior of these correlators. 

To obtain an action suitable for describing phase (I), we 
apply the duality procedure outlined in the Appendix to the 
Ji variables in our initial action. We obtain the following 
reformulation;^ 



S[Qx,J 2 ] = \J2 
k 



\2irQt(k) + 0(k)J 2 (k)\ 



*'i(fc)IM 2 
(2-) 2 



vi(k)\fk\ 




vi(k)\f k \ 2> 



y(-k)-J 2 (k) 



(14) 



where fk tfl = 1 — e 4 ^, as defined in the Appendix. The ac- 
tion is written in terms of J 2 variables, and Q\ variables that 
are dual "vortex" variables to J\. In the above action, and 
from now on, we consider the case of general intra-species 
interactions v\(k) and v 2 {k), though in the preceding section 
we considered specific short-range interactions v\(k) = 1/ti 
and V2(k) = l/t 2 . We also consider a more general, k- 
dependent statistical coefficient 9(k). Throughout this work, 
we will assume that v\ (k),v<2(k) and 9(k) are real and satisfy 
v a (k) = v a (—k) and 9(k) = 6(—k). We can see in the above 
action that for large t\ and small t 2 , both Qi and J 2 have a 



large energy cost, and therefore both are gapped. We expect 
this; since the J\ variables are condensed, the variables dual 
to them should be gapped. Naturally, we can obtain a refor- 
mulation for phase (II) by applying the same steps to the J 2 
variables. 



To get an action suitable for describing phase (IV), we ap- 
ply the duality procedure to the J2 variables in Eq. (fl4l . This 
gives us the following action, expressed in terms of "vortex" 
variables Qi and Q 2 that are dual to the J\ and J 2 variables: 



J 



(2tt) 2 v 2 (k)\Q 1 (k)\ 2 +v 1 (k)\Q 2 (k)\ 



|/fe| 2 ^l(fc>2(fc)+0(fc) 2 



r 



{2nfe(k) Qi(-fc) -PQ2(k) 
|/ fc | 2 «i(fc)« 2 (fe)+0(fc) 2 



(15) 



Here pq 2 is an auxiliary "gauge field" encoding the flux of 
Q 2 , and is defined such that Q 2 = V x pq 2 (because of the 
constraints on Q1.2, the action is independent of the choice of 
Pq 2 ). Unlike the analysis of phase (I) in the Qi, J 2 variables, 
it is not clear that a phase with gapped Q\, Q 2 exists. If we 
define i>i/2,dual sucn m at 



Ul/2,dual(*) 



(27T) 2 «2/l(fc) 



\fk\ 2 vi{k)v 2 {k) + e{kf 



(16) 



we see that fi/2,duai cannot both be arbitrarily large, and so 
the interactions may not be large enough to gap out both Qi 
and Q 2 . Considering comparable V\ ~ v 2 , the dual interac- 
tions are largest for intermediate V\ and v 2 and their magni- 
tude increases with decreasing 9. Whether a phase with both 



species gapped exists needs to be determined numerically. We 
have found that in the current model with = 2ir/3, phase 
(IV) is the phase where Q\ and Q 2 are gapped. In contrast, in 
the 6 = it model^ such a phase did not exist, and instead we 
expect either a critical state or phase separation.i2i2li 

We also note that the transition from phase (IV) to phase (II) 
is a transition where the Qi variables are going from gapped to 
condensed, while the Q 2 variables remain gapped. Therefore, 
if we could study correlators such as (Qiu(k)Qiu(— k)), we 
would expect them to behave in the well-understood manner 
of one species of loop condensing, qualitatively similar to the 
behavior of the J 2 variables in the (O)-(II) transition. We now 
invoke a useful relation derived by introducing external probe 
fields as explained above: 
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C 22 J(k = (0,0, k z )) 



vm\h 



\h\ 2 vi{k)v 2 (k)+9(kf 



(17) 



(2^) 2 


^(fc) 2 (Oi,(fc)gi y (-fc))-| 


f k \ 2 v 1 (k) 2 (Q 2x (k)Q 2x ( 


-*)> + 


Asm{k z /2)v l {k)9{k){Q ly {k)Q 2x (~k)) 




)fk\ 2 v 1 (k)v 2 (k) + 


e(k) 2 


2 



Q2 is gapped everywhere near the (IV)-(II) phase transi- 
tion, which implies that the excitations of Q 2 are small 
loops, and we can show that (Q 2fl (k)Q 2 ^(—k)) ~ k 2 for 
small k in the region of the transition. We also expect that 
(Qiy(k)Q 2x (—k)) ~ fc 3 in phase (IV) and ~ k in phase (II). 
Taking the limit of small fc, we see that most of the terms are 
of order fc 2 or smaller, and we are left with 

C l» = S<<M*)Qiv(-*)> + °( fc2 )- ( 18 > 
This explains why C 22 (fc m j n ) ~ 1/L 2 in phase (IV) and is 



constant in phase (II), even though J 2 is condensed in both 
phases. It allows us to use C 22 to study the condensation of 
the Qi variables, which is what we showed in Fig. [3] Equa- 
tion ( fT8l is also valid at the transition between phase (IV) and 
phase (III), shown in Fig.[lO| We can see from this figure that 
the Q variables seem to be undergoing a continuous transition; 
we will discuss this further in Sec. IV B| 

We can also establish the behavior of C xy (k) in phase (IV) 
by invoking another relation, again derived by differentiating 
the partition function with respect to external probe fields : 



2sin(ft*/2)0(fc) 
\fk\ 2 vi(k)v 2 (k) + 8(ky 



(27r) 2 2sin(fc z /2)6»(fc) 



(2^) 2 


"|/fc|^l(%2(fc)-0(fc) 2 " 






\ 2 V!(k)v 2 {k) + <9(fc) 2 


2 



\fk\ 2 V!(k)v 2 (k) + 9(k) 2 

{Qx y (k)Q 2x (-k)}. 



[vi(k)(Q 2x (k)Q 2x (-k)) + v 2 (k)(Q ly (k)Q ly (-k))} 



(19) 



We can see that for gapped Qi. 2 and in the small fc limit, 
C x 2 indeed approaches —k/0, as we argued from a schematic 

treatment of the Ji, J 2 condensate in Eq. ( TT2l . Note that we 
can use Eqs. ( TTTb and (19[ to express the Q correlators in 
terms of the J correlators. We have done this, and plots of 
the data (not shown), confirm the condensation of Qi across 
the (IV)-(II) and (IV)-(III) transitions, as expected from the 



above analysis. We chose to express all data in terms of corre- 
lators in the J variables so that Section|III]could be understood 
without any knowledge of the various reformulations. 

We can also give precise meaning to the treatment in 
Eq. dT2b . From the Appendix, an intermediate step in the du- 
ality procedure going from Ji, J 2 to Qi, Q 2 is: 



S[dji 7 dj 2 ,Qi,Q 2 ] = 



1 ^ [V x aji(R)] 2 



E 

R 



(2^) 2 ii 



v [V x a j2 (r)} 2 
4- (2n)H 2 



(2tt): 



5^[V x aj^R)} ■ a j2 (r) 



+iJ2Qi(R)- Sji (R) + i E (r) • ct j2 (r) . 

_R r 



(20) 



r 



Gaussian integration over the a variables gives Eq. (TT3T >. Chem-Simons interactions coupled to integer-valued currents 
Equation (|20]) is an action for two gauge fields with mutual Q. When Qi and Q 2 are gapped, we can formally integrate 
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them out and obtain the low-energy field theory description in 
Eq. Q3. 

We now consider a reformulation appropriate for the de- 
scription of phase (III). Our crude intuition is that the J\ and 
J 2 loops will indeed condense strongly but only in multiples 
of n, where 9 = — , in order to avoid the statistical inter- 

action. To proceed more accurately, we start with S[Qi, J 2 ], 
Eq. ( TBI ), and notice that for small vi and v 2 the combination 



M 2 (R) = J 2 {R) + nQ 1 (R) 



(21) 



wants to be gapped while Qi wants to be "condensed", hence 
J 2 wants to be "condensed" in multiples of n. More precisely, 
in the partition sum, we can change the summation variables 



J 



from integer- valued currents Qi and J2 to integer- valued cur- 
rents Qi and M 2 , with the action 



(2n) 2 \M 2 (k)f 



(22) 

n 2 vi(k)\f k \ 2 
i^i; 2 (fc)|M 2 (fc)-nQi(fc)| 2 . (23) 



We can now consider a phase with M 2 gapped and Qi con- 
densed appropriate for small v\ and v 2 . The precise meaning 
of the Qi condensation is again obtained by going from Qi 
to dual variables Mi using the formal prescription in the Ap- 
pendix. The result is 



S(e=W")[M 1) M 2 ]=X]-^ ! 2 



» 2 \U 



Wi{k)\ 2 t |M 2 (fc)p 



2v 2 {k) 2v 1 (k) 



;V— M 1 (-k)-p M 2(k) 

k 



(24) 



r 



where M 2 = V x f>M 2 . Note that if we were to dualize Qi 
in S[Qi, J 2 ], we would of course obtain back S[Ji, J 2 ] (up 
to sign of Ji), while the duality procedure after the change 
of variables in Eq. (fJTJ gives a different reformulation since 
here we dualize Qi while keeping M 2 as an independent cur- 
rent. Labels 1 and 2 on Mi and M 2 are somewhat arbitrary, 
as we mixed the original species 1 and 2, e.g., when defining 
M 2 in Eq. (fJTJ. We can think about a phase with gapped M 2 
as having binding of n original J 2 currents to one anti-vortex 
in Ji, so that for J 2 = n, Qi = — 1, we have M 2 = 0. In 
phase (III) it is such (J 2 = n, Q\ = —1) molecules that are 
condensed. These bound states are illustrated in Fig.QT] This 
is the more accurate description of phase (III), made precise 
by the reformulation Eq. (l24l with gapped Mi. 2- For n = 2, 
this is also the precise description of phase (III) in Fig. 1 of 
the 6 = 7r statistics model studied previously.^ We can treat 
a small loop of M 2 as an excitation out of such a bound state. 
The symmetric structure of the above action suggests a simi- 
lar interpretation of Mi even though this field was introduced 
differently. 

Since the M variables are gapped in phase (III), we know 
from small loop arguments that (M afJi (k)M a ^(—k)) ~ k 2 and 
(Mi AI (fc)A/2i/(— k)) ~ fc 3 . We can also derive the following 
exact relations at k = (0, 0, k z ): 



c: 



1 



(27T) 5 



v 2 {k) n 2 \f k \ 2 v 2 (k) 

(27T) 2 



■{Mi x {k)Mi x (-k)\25) 



C 12 = 

xy e - 1 2 



^-(Mi x (k)M 2y {-k)) (26) 



n 2 vi{k)v 2 {k)\f k \ 
These imply that C 22 (k) -(constant) and C 12 {k) — k. Note 



that unlike in phase (IV), C 12 (/c m ; n ) ■ L has a non-universal 
value in phase (III), as can be seen in Fig. [5] 




FIG. 1 1 : Illustration of the "molecules" that are condensed in phase 
(III). Each molecule contains three charges (stars), and one anti- 
vortex (circles). Black objects on white backgrounds are of species 
1, and white objects on black backgrounds are of species 2. 



V. FIELD THEORIES FOR THE PHASES AND PHASE 
TRANSITIONS 

A. Phase (0) and the (0)->(IV) transition 

Equation ( fT2b is a continuum field theory useful for describ- 
ing phase (IV). In this phase, the J variables are condensed 
which allowed us to use real-valued variables j = 



Vxa, 



2tt ' 

In phase (0), the Q variables are condensed, and we can re- 
place them with real-valued variables q. We can then write a 
field theory in terms of real-valued gauge fields a q such that 

' Performing this procedure on Eq. (fT31 gives: 



q = 



2 71 



10 



S[d q i,a q2 , Ji, J 2 ] 



k 



k L 



V2(k)\[V x a ql ](k)\ 2 
fk\ 2 v 1 {k)v 2 {k)+9{k) 2 
6(h) 



«i(fc)|[v x d, 2 ](fc)| 5 



\fk\ 2 vi(k)v 2 (k) + 9(kf 



(27) 



IM 2 «i(£> 2 (fc) + W 



[V x d q x](-k) ■ a q2 {k) |/M~ fc ) ' ®qi( k ) + M-k) ' « g 2(fc) 



r 



which can be viewed as an intermediate step in the (exact) 
duality map from the variables Q1.2 to Ji. 2 , cf. the Appendix. 

I 



We can now take the long-wavelength limit and write a 
schematic action in real space: 



Soft [d q l, d q2l ^,12} = I d 

d 3 r 



(V x d ql f (V x d q2 f i _ - _ 

Q a i 1 ' ( V X a i 2 > 



2t 2 2 2t 1 9 2 

71 1 (V - iaqOfjil 2 + 72 1 (V - id q2 )^j2\ 2 + mi|*ji| 2 + m 2 |*,/ 2 | 2 + (quartic terms) 

I 



(28) 



where we used continuum complex- valued fields ^fji, v Pj 2 
to represent the matter that was represented on the lattice by 
the integer- valued currents J x , J 2 , and we did not write the 
quartic terms explicitly. This is the action for two gauge fields 
with mutual Chern-Simons interactions, and two matter fields, 
minimally coupled to the gauge fields. Here 71, 72, mi, m 2 
are some effective parameters; along the self-dual line we 
have 71 = 7 2 and mi = m 2 . For gapped ^fji, v I'j 2 , we 
can integrate these out and obtain a long-wavelength descrip- 
tion of (0) in terms of the a q variables. Condensation of 
^>j2 leads to phase (IV). We therefore propose Eq. (f28) as the 
field theory describing the transition at the lower multicritical 
point. As discussed in Sec. [TTIJ our results on the nature of the 
(IV)-(O) transition are still conflicting, but we hope that they 
will stimulate further numerical and analytica l 10 i 47 ~ — studies. 



B. Phase (III) and the (III)^(IV) transition 

To get another perspective on phase (III), we first interpret 
the coefficient on the last term of Eq. (fT3T > as a statistical in- 
teraction for the Q variables, given by 



7dual 



(k) 



-(2n) 2 9{k) 



\fk\ 2 v 1 (k)v 2 (k)+9 2 (k) 



(29) 



We can shift this coefficient by 2itn, for integer n, without 
changing the Boltzmann weight e~ s . This gives us the fol- 
lowing equation for the new statistical angle: 



-'dual, shifted 



(k) 



(27T) 2 \f k \ 2 V 1 (k)v2(k) 



9 \f k \ 2 Vl (k)v 2 (k) 



(30) 



where we have used 9 = 2tt/ii. Performing formal duality 
on Eq. ( fT5l ) using the new statistical interaction 0duai,sMfted 
gives precisely Eq. (|24| |. [The non-commutation of the du- 
ality and shift of 9 by multiple of 2tt is well known in the 
literatur e) 29 ' 32-34,40 and is known to correspond to possibility 
of "oblique confinement". Here we explicitly see this relation 
by identifying precise bindings of objects in phase (III) as de- 
scribed in Sec. [IV]] This allows us to interpret the variables M 
that are gapped in phase (III) as being dual to the Q variables 
after the shift. The precise meaning of the duality is as in the 
Appendix, and can be viewed as replacing the integer-valued 
Q variables by real-valued variables q, while maintaining the 
information about integer-valuedness in terms of new integer- 
valued M. We can write q = — l^- 3 , and obtain the following 
action: 



S ^M0 ql J q2 ,M 1 ,M 2 } = lj2 



vi{k)v 2 (k) 



v 2 (k)\[W x (3 ql ](k)\ 2 Vl {k)\[V x q2 ](k)\' 



\fk\ 2 vi(k)v 2 (k) + 2 |/ fc | 2 «i(fcMfc) + 9 2 



+ ^ 0[\f k \ 2 Vl (k)v 2 (k) 



9 2 



[V x V x p ql ]{-k) • [V x (3 q2 ](k) + 1^ [Mx(-k) ■ gl (k) 



(31) 

M 2 {-k) ■ p q2 (k) . 
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Compared to Eq. ( f2Tb . we have used different labels for the 
gauge fields even though the first terms are the same, to em- 
phasize that the coarse-graining procedure will have a dif- 
ferent meaning after the shift in #duai- Again, we can cast 
this action into a more familiar form by returning to real 



space while taking the long-wavelength limit and replacing 
the current-loop representation with complex matter fields 
^mi and *a/2, minimally coupled to (3 q i and (3 q2 : 



-,(e=27r/ra)r 
'off I 



S a(r \Pql,Pq2 



d 3 T 



M2 



(Vx/3 9l ) 2 (Vxfef 



+ S77(VxW-(VxVx^) 



2t 2 6 2 2t 1 9 2 BH^ 

71 1 (V - z/V^a/iI 2 + 72 1 (V - ^2)^/2 1 2 + mi |*Afi | 2 + m 2 \^ M2 \ 2 + (quartic terms) 



(32) 



In phase (III) the $mi. fields are gapped and we 

can integrate them out to obtain a description in terms of 
two gauge fields j3 q i and j3 q 2, with a higher-order mutual 
Chern-Simons term^^ The latter is irrelevant compared to 
the Maxwell terms and does not gap the gauge fields. We 
thus have two gapless modes. The transition from phase (III) 
to phase (IV) is a condensation of the ^mi and \Pm2 vari- 
ables, and along the self-dual line the two fields condense 
simultaneously. We conjecture that the higher-order mutual 
Chern-Simons term is irrelevant at this transition as well, 
and therefore the transition is two decoupled inverted XY 
transitions! 37 ' 53 ' 54 Returning to Qi, Q 2 variables, we conjec- 
ture that the transition from phase (IV) to phase (III) is two 
decoupled XY transitions where the mutual statistical interac- 
tion of Qi and Q 2 , given by 8 dua i sh iftcd in Eq. d30j is irrel- 
evant at long wavelengths. This interpretation, together with 
Eq. ([T8]l, explains why we are able to use finite-size arguments 
on the data in Fig.[lO]to study the properties of the (IV)-(III) 
phase transition and conclude that it is continuous. If this in- 
terpretation is correct, it also means that the phases (IV), (III), 
(I), and (II) all meet at a single multicritical point. 

We remark that while the lattice actions Eqs. ((27) and ( |3D 
are mathematically equivalent and contain all phases in Fig. [2] 
the continuum field theories Eqs. (l28l l and (l32l are distinct and 
apply only near the corresponding multi-critical points. 



VI. IRREDUCIBLE RESPONSES 

The current-current correlators C c ^ v represent the response 
of the current J afJi to an externally applied field A^f. In sys- 
tems with long-range interactions it is useful to study "ir- 
reducible responses" C^ ,lrred , which are the responses of 
the currents to the total field A tot , made up of both A ext 
and an internal gauge field induced by the other currents in 
the system^— In our model, the statistical interaction is 
the long-range interaction, and it acts between different loop 
species in perpendicular current directions. In this section we 
will derive the appropriate expressions for the irreducible re- 
sponses, and show their behavior in our system. 

If we apply external fields coupled to both species of loops, 



as in Eq. ( fT3l ), then by the Kubo formula the response of the 
current variables is given by: 



{Ja li {k)) = -iY,Cf v {k)AZ\k). 



(33) 



For concreteness, we will assume that k is in the z direction, 
k = (0, 0, k z ), and this implies that = if \i or v are 
in the z direction. As discussed in Sec. [II] the lattice mirror 
symmetries of our action mean that the only correlators which 
are non-zero in Eq. (|33l are C"" and C^ 2 with /i ^ v. This 
implies that for a gauge field in one direction, we need only 
to consider its effects on two of the six possible J afi ; for con- 
creteness in this work we will consider J\ x and J 2y . This 
allows us to write Eq. ( |33l in the following way: 



<J> = 



(Ji x ) 



C = 



J) = -iCA cxt 



(34) 



Til 

J XX 

-Cll C? 



■'■// 



mi 



A CXt 

A lx 

A CXt 

A 2y 



where we have used the fact that C, 2 ^. = — Cj, 2 , which we 
can also deduce from the mirror symmetries of our model. To 
characterize the response of (J) to the total field, we write 



-iC 



irrcd a tot 



(35) 



with C llrcd and A tot defined similarly to the quantities in 
Eq. ( f34b . Here A^ is the total field, and is identified as 



A, 



(36) 



where the gauge fields a qafl are precisely those in Eq. ( |27i i 
mediating the J\ and J 2 interactions. We can calculate the 
expectation values (a q ) in the presence of A cxt by analyzing 
Gaussian integrals in Eq. ( |27] i: thus, for any fixed J we obtain 
the following: 



(a) = iVJ 

(ai x ) 
(a 2y ) 



(37) 



V 



vi(k) 
-g(fc) 



2sin(fcj/2) 
V 2 (k) 



Inserting this into Eq 

^tot 



2sin(fc^/2) 

and using Eq. (l34l we get 

(l-VC)A 0Xt , (38) 
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and comparing Eqs. ( l34t and ( l35l l gives our final expression 
for the irreducible responses: 



C irr 



C(l- VC) 



(39) 



We can use the irreducible responses to determine the con- 
ductivities of the system, through the relation 



(40) 



We have plotted the diagonal and off-diagonal conductivities 
along the self-dual line in Figs. [T2l and [T3l As in Ref. H3, we 
can use these conductivities to detect condensation in systems 
with long-range interactions. We can see that a 22 diverges 
with the system size and thus detects the condensation of J\, 
J 2 in phase (IV), while we recall from Fig.|4]that the correlator 
C 22 did not. The diagonal conductivity has a crossing at t\ 
t 2 ~ 0.338, which is tentatively the location of the transition 
between phases (0) and (IV). 

It is interesting to note that in phase (III), <r 22 = while 
al 2 approaches a universal value of 1/6 = 3/(271"). We can 
loosely interpret this if we recall that phase (III) is a conden- 
sate of composite objects containing n particles of one type 
bound to one antivortex of the other type [see Fig. Q~TJ and 
Eq. d2"Tli1. For example, consider a situation where we have 
a J\ charge current flowing in the x-direction. This current 
can be carried by the condensate of the bound states, in which 
case there is also a Q 2 current in the x-direction, given by 
(Q2x) = —(Jix)/n- In the absence of any other currents, 
we get A\°* — 0. Furthermore, we can think of the Q 2 
variables as magnetic fluxes for the J 2 charges. Therefore, 
by Faraday's law there is an electric field (acting on the J 2 
charges) induced perpendicular to the direction of Q 2 , and 
we get -ik z A 2 °y = -2tt(Q 2x ) = ^Jix- This is exactly 
what we would expect from the conductivity that we derived, 
A tot = i (c irred )- 1 (J). 

We can consider the responses for the dual Q\, Q 2 variables 
as well. Focusing on a pair Q\ y , Q 2x , we define 



^dual 



(Qi v (k)Q ly (-k)) (Q ly (k)Q 2x (-k)) 
-{Qiy{k)Q 2x (-k)) (Q 2x (k)Q 2x {-k)) 



where we used (Q2x(k)Qi y (—k)) = -(Qi y (k)Q 2x (-k)). 
The interaction matrix for the specific ordering of the cartesian 
components is 



V 



dual 



^l,dual(fc) 
8du„i(fc) 



-9d„al(fc) 

2sin(fc i /2) 



2 S in(t/2) «2,dual(fc) 



\fk\ 2 



(y- 1 ) 1 , (4i) 



where the last relation was derived by using Eqs. ([Tol and 
j29l , and the superscript "T" denotes the matrix transpose. 
The dual and direct responses satisfy the relation 



vc + c 



.al^dual — lj 



(42) 



which we can check by using Eqs. ( fTTI ) and ( fT9l . Re- 
lation ( f42b is similar to the relation satisfied in the one- 
component casei^t^ We can also verify that the irreducible 




= FIG. 12: The conductivity a , along the self-dual line, obtained 
from the raw data in Figs. |4] and [5] Vertical lines mark the phase 
boundaries. We can see that a 22 diverges in phase (IV), but is zero 
in phase (III). 
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FIG. 13: The transverse conductivity a^. 2 , along the self-dual line, at 
the boundary between phases (IV) and (III). A vertical line marks the 
phase boundary. In phase (IV) a 12 y approaches a non-universal value, 
while in phase (III) it approaches the universal value 3/2n, which is 
indicated by the horizontal line in the figure. Calculations of aiy in 
phase (IV) involve cancellations of similar quantities, which greatly 
increases the noise in this region. The origin of the cancellation is in 
the quantization of the C 12 values in this phase, as seen in Fig. 1 101 
Therefore we have chosen not to display our data in phase (IV). 



conductivities satisfy 



"""dual - 



(27T)S 



(43) 



which is similar to the relation that conductivities obey when 
there is only one species of loop^r— 
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VII. DISCUSSION 

It is instructive to compare these results to those of our ear- 
lier study at 9 = ttJ£- The Boltzmann weight e~ s of the 9 = tt 
model is invariant under (Jx, J2) — > (—Ji, J2), while this is 
not satisfied in the present model. We can see that the cor- 
relation between different currents, C 12 , changes sign under 
this operation, and therefore must be zero in the 9 = tt model. 
This explains why that model did not contain the phase (IV) 
that we have seen in the present study. The location of phase 
(III) in the two models is also quantitatively different, since 
the loops must condense in different strengths to avoid the 
statistical interaction, and this happens at different values of 
t. Phase (III) itself is qualitatively similar in the two models 
(except for the charge multiplicity in the bound states), and is 
detected by (7 12 > lrred i n the 9 = tt model despite the fact that 
C 12 is strictly zero. 

From these studies, we can anticipate the behavior of the 
model with short-range interactions at general statistical an- 
gle 9 7^ tt. We expect that the phase diagram will be similar 
to the one in Fig. [2] except that the phase (III) will feature 
condensation of more complex composites^— and will oc- 
cur at different values of t. An open question in the present 
work is the nature of the lower multi-critical point, where our 
results are conflicting between first and second-order scenar- 
ios. It would be interesting to explore this phase transition 
in short-ranged models for other values of 9 numerically and 
analytically. 

It is also interesting to explore behavior for more general 
interactions, particularly for self-dual models with Ji <H- J2 
interchange symmetry. For the model with short-range inter- 
actions, we have seen that the statistical interaction qualita- 
tively changes the nature of the phases and phase transitions. 
On the other hand, for loops with long-ranged interactions de- 
caying as 1/r in real space (behaving as 1/fc 2 for small k 
in momentum space), we expect that the statistical interac- 
tions are less important, since here the density fluctuations are 
very strongly suppressed and the mutual statistics phases are 
fluctuating lessi 5 ^ 5 ^ In fact, starting with the original model 
Eq. ([l} with short-ranged interactions at 9 = 2n/n, our re- 
formulation in terms of Mi and M% variables in Eq. (l24l can 
be viewed as precisely such a new model with long-ranged 
interactions and 9mi.M2 = 27r/n, so the present numerical 
study already provides information about such a model with 
$mi,M2 = 27r/3. In the absence of the statistical interac- 
tions, loops with long-ranged interactions would condense via 
independent one-component Higgs transitions (inverted XY 
transitions) ^ 7,53 ' 54 From our discussion in Sec.|V] we conjec- 
ture that this remains true also in the presence of the statistical 
interactions with 9 ^ tt, i.e., they are irrelevant at the phase 
transition in the long-ranged case. 

An interesting case is obtained for marginally long-ranged 
interactions decaying as 1/r 2 in real space (behaving as l/|fc 
for small k in momentum space).^ 5 ^ In a recent paper Id7I1 . 
we studied condensation of single species with such marginal 
interactions and found second-order transitions with continu- 
ously varying critical properties that depend on the coupling 



of the long-range interaction. We would like to study conden- 
sation for two species with mutual statistics and ask whether 
the transitions remain continuous for 9^0 and explore the 
critical properties (which will likely vary with 9). We can 
construct a lattice model where we know the phase boundaries 
exactly from duality considerations^— and can focus on such 
studies precisely at the transitions. An interesting question is 
what happens for 9 = tt in such models with marginally long- 
ranged interactions, whether we find a critical loop state or 
phase separation. The latter happened in a specific model with 
short-ranged interactions that we studied in Ref . HH while we 
would like to explore if a critical state can be obtained for 
modified interactions. 

For broader outlook, our system is an example where cer- 
tain reformulations allow direct study of particles with mutual 
statistics. It would be interesting to look for other cases where 
such reformulations may be possible. Systems with more 
complex anyons could be interesting ^ 6,7,13,15- — , and such com- 
bined numerical and analytical studies could bring insights 
about broader phase diagrams and phase transitions involving 
topological phases. Furthermore, the present two-loop sys- 
tem can be viewed as an example of more general actions 
with topological terms. In fact, as discussed in Ref. Qjl the 
two-loop model with 9 = tt statistical interaction is equiv- 
alent to an anisotropic 0(4) sigma model with a topological 
9 = tt term; our loop models can be viewed as providing 
precise lattice realization of this topological field theor y 18 ' 28 
and show that it is important to examine different phases 
such a theory may have. Inspired by our two-loop systems, 
it would be interesting to study precise lattice (discretized 
space-time) formulations of other topological field theories of 
current interest^^ r 16 ' 25,27,60- — also in other space-time dimen- 
sionalities, and ask if they may also allow sign-free reformu- 
lations and hence unbiased numerical studies. 



Acknowledgments 

We are thankful to A. Vishwanath, M. P, A. Fisher, 
T. Senthil, R. Kaul, G. Murthy, J. Moore, N. Read, A. Shapere, 
and W. Witzak-Krempa for stimulating discussions. This 
research is supported by the National Science Foundation 
through grant DMR-0907145; by the Caltech Institute of 
Quantum Information and Matter, an NSF Physics Frontiers 
Center with support of the Gordon and Betty Moore Founda- 
tion; and by the XSEDE computational initiative grant TG- 
DMR 110052. 



Appendix A: Formal duality procedure 

This appendix summarizes our duality procedure for one 
loop species i 11,35- — The original degrees of freedom are con- 
served integer-valued currents J(r) residing on links of a sim- 
ple 3D cubic lattice; V • J(r) = for any r. To be precise, we 
use periodic boundary conditions and also require vanishing 
total current, J tot = ^ r J(r) = 0. We define duality map- 
ping as an exact rewriting of the partition sum in terms of new 
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integer-valued currents Q(R) residing on links of a dual lat- 
tice and also satisfying V • Q(R) = for any R and Q tot = 0: 



E 1 



-SduallQ] 



(Al) 



-S dual [Q=Vxp] 



(A2) 



In the first line, the primes on the sums signify the above con- 
straints on the currents J and Q respectively. In the second 
line, the prime on the real-valued integration measure signi- 
fies corresponding linear constraints realized with Dirac delta 
functions, Il r: ^o5[V ■ j(r) = 0] and <5(jtot = 0). For any 
configuration Q satisfying the above constraints, we can find 
p{r) such that Q = V x p, and the constraints on j guarantee 
that the right-hand-side of the last equation does not depend 
on the choice of p. 

Equations ( IAU - (IA2b provide a precise way to go from 
integer-valued sums with constrained J to real-valued inte- 
grals with constrained j, which is achieved with the help of 
new integer-valued constrained fields Q. A formal demon- 
stration can be sketched, e.g., as follows: We first im- 
plement the constraints on J using conjugate 27r-periodic 
phase variables. We then replace sums over integer-valued 
Jfi{r) with integrals over real-valued j^(r) containing a fac- 



tor£ 



-«V(r)27rp M (r) f Qr eacn j m j, ^ e g roU p con . 



figurations p(r) into classes specified by Q = V x p and use 
summation over members in each class to effectively extend 
the integrations over phase variables to the full real line. The 
latter integrals finally lead to the delta function constraints on 
the real-valued fields j defining the measure [Dj]'. In the pro- 
cess, we see that Q can be interpreted as vortex lines in the 
phase variables conjugate to J. 

An immediate important application is to the case with 

= ^E l, ( fc )l / ( fc )| 2 + l E J(-k)-A^(k), (A3) 



where we have also coupled the original currents to an exter- 
nal probe gauge field A cxt . The integration over j in Eq. iA2i 
is Gaussian and readily gives basic averages 

iUk)jAk'))o= S -^(s^-^^) , (A4) 



v(k) 



where fk 4L = 1 — e >* . We then obtain 



\fk 



- 1 ^ [2nQ(-k) + B(-k)} ■ [2nQ(k) + B(k)} 
1 k v(k)\fk\ 2 

(A5) 

where B = V x A cxt . The relation between Eq. jA3i and 
Eq. (|A5l l is what we call "duality map" in the main text. 
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